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Abstract. We present a detailed theoretical and experimental study of Feshbach resonances in the 6 Li- 40 K 
mixture. Particular attention is given to the inelastic scattering properties, which have not been considered 
before. As an important example, we thoroughly investigate both elastic and inelastic scattering properties 
of a resonance that occurs near 155 G. Our theoretical predictions based on a coupled channels calculation 
are found in excellent agreement with the experimental results. We also present theoretical results on the 
molecular state that underlies the 155 G resonance, in particular concerning its lifetime against spontaneous 
dissociation. We then present a survey of resonances in the system, fully characterizing the corresponding 
elastic and inelastic scattering properties. This provides the essential information to identify optimum 
resonances for applications relying on interaction control in this Fermi-Fermi mixture. 
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1 Introduction 

A new frontier in the research field of strongly interact- 
ing Fermi gases [1,2] has been approached by the recent 
realizations of ultracold Fermi-Fermi mixtures of 6 Li and 
40 K [3-9]. Degenerate Fermi- Fermi mixtures represent a 
starting point to experimentally explore a rich variety 
of intriguing phenomena, such as many-body quantum 
phases of fcrmionic mixtures [10-24] and few-body quan- 
tum states [25-28]. 

The possibility to precisely tune the interspecies inter- 

\ action via Feshbach resonances [29] is an important pre- 
requisite for many experiments. This has motivated the- 
oretical and experimental work on Feshbach resonances 
in the 6 Li- 40 K mixture [4,7,30]. It turned out that all 

: resonances for s-wave scattering in this system are quite 
narrow, the broadest ones exhibiting a width of <2G 1 , 
and their character is closed-channel dominated [29] . This 

' causes both practical and fundamental limitations for ex- 
perimental applications. Interaction control is practically 
limited by magnetic field uncertainties and, on more fun- 
damental grounds, the universal range [29] near the center 
of the resonance is quite narrow. 

Our work is motivated by identifying the Feshbach res- 
onances in the 6 Li- 40 K system that are best suited for re- 
alizing Fermi-Fermi mixtures in the strongly interacting 
regime. In a previous study [7], Tiecke et al. approached 



Here, G= 0.1 mT. 



this question by calculating the widths of the different res- 
onances [31], and they studied elastic scattering for one of 
the widest resonances in the system. Another important 
criterion is stability against inelastic two-body decay. For 
the 6 Li- 40 K mixture, inelastic spin-exchange collisions do 
not occur when at least one of the species is in its low- 
est spin state [4]. When one of the species is in a higher 
state, decay is energetically possible, but rather weak as 
it requires spin-dipolc coupling to outgoing higher partial 
waves. The wider resonances in the 6 Li- 40 K system are 
found in higher spin channels. This raises the important 
issue of possible inelastic two-body losses. The question 
of which is the optimum resonance for a particular appli- 
cation can only be answered if both width and decay are 
considered. 

In this Article, we present a detailed study of Feshbach 
resonances in the 6 Li- 40 K system, characterizing their in- 
fluence on both clastic and inelastic scattering properties. 
In Sec. 2 we briefly review a general formalism to describe 
decaying resonances. In Sec. 3 we present a case study of 
a particularly interesting resonance. In theory and exper- 
iment, we investigate its clastic and inelastic scattering 
properties and the properties of the underlying molecular 
state. In Sec. 4 we present a survey of all resonances, sum- 
marizing their essential properties. In Sec. 5 we conclude 
by discussing the consequences of our insights for ongo- 
ing experiments towards strongly interacting Fermi-Fermi 
mixtures. 
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2 Feshbach resonances with decay 

A Feshbach resonance results from the coupling of a col- 
liding atom pair to a near-degenerate bound state. If this 
molecular state can decay into open channels other than 
that in which the colliding pair is initially prepared, the 
situation is referred to as a decaying resonance [29] . A for- 
malism for describing such resonances has been developed 
for optical coupling [32,33], and has also been applied to 
the magnetically tunable case [34]. A well known exam- 
ple of a decaying resonance exists in the collision of two 
85 Rb atoms [35], for which molecular lifetimes have been 
studied [36,37]. 

The scattering properties in the zero-energy limit can 
be expressed through a complex s-wave scattering length, 
a = a — ib, where a and b are real. The relation of these 
two parameters to the two experimentally relevant quan- 
tities, the elastic scattering cross section a and the loss 
rate coefficient for inelastic decay, is for non-identical 
particles given by 



4ir(a 2 + b 2 ) 



and 



K' 



- 2k b 



(1) 



(2) 



Here, h is Planck's constant and /i is the reduced mass. 
The complex scattering length can be parametrized by 



a(B) = a hg 



b(B) = 2a res 



7b (-B — Bo) 



(B - So) 2 
(7b/2) 2 
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(4) 



Here, B is the magnetic field strength, the resonance 
occurs at B = B$, and o,b g is the background scatter- 
ing length. The decay of the "bare" molecular state that 
causes the resonance is characterized by a rate 7 [37], 
which we conveniently express in magnetic field units, 
7b = fvy/Sfi, where 5[i is the difference in magnetic mo- 
ment between the entrance channel and the bare molecu- 
lar state. The resonance length parameter a ros is related 
to the resonance width A by 



7b 



(5) 



and gives the range ab g ± a r os within which the real part 
of the scattering length can vary, thus providing an indi- 
cation of the possible control. A common figure of merit 
for the coherent control of an ultracold gas is the ratio 
a/b. For \B — Bq\ ^> 73, and a change in scattering length 
much larger than ab g , this can be shown from Eqs. (3) and 
(4) to be 



, (B-Bq) 

7b 



(6) 



A larger a res therefore gives better coherent control and 
lower losses for a given change in scattering length. Com- 
bining Eqs. (2) and (6) gives a simple expression for the 
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Fig. 1. Zeeman sub- levels in the electronic ground state 
of 6 Li and 40 K, giving the total angular momentum / and 
its projection mf along the quantization axis. We label 
Zeeman states alphabetically in order of increasing energy, 
as shown. 



loss rate coefficient 



K 2 « — 



(7) 



which shows a general a 2 -scaling of two-body loss near a 
decaying Feshbach resonance. 



3 Case study of the 155 G resonance 

In this Section, we present a thorough study of the 155 G 
Feshbach resonance, which serves as our main tool for in- 
teraction tuning in strongly interacting Fermi-Fermi mix- 
tures. It was first observed in Ref. [4] and used for molecule 
formation in Ref. [5]. We first (Sec. 3.1) present theoretical 
predictions for the elastic and inelastic scattering proper- 
ties near this resonance based on coupled channels cal- 
culations. We then (Sec. 3.2) present our corresponding 
experimental results, providing a full confirmation of the 
expected resonance properties. We finally (Sec. 3.3) dis- 
cuss the properties of the molecular state that causes the 
resonance. 

Figure 1 shows different magnetic and hyperfinc sub- 
levels of the electronic ground states of 6 Li and 40 K. Here 
we follow the notation of Ref. [29] and label the sub-states 
alphabetically in increasing order of energy. The 155 G 
resonance occurs in the ac channel, i.e. for a 6 Li atom in 
state a colliding with a K atom in state c. 



3.1 Scattering properties: Theory 

We have carried out coupled channels studies [38] of the 
scattering properties in the ac s— wave channel. The po- 
tentials used were taken from Ref. [30] and, to make this 
paper self-contained, we have summarized the important 
parameters in Table 1. The resonance is created by spin- 
exchange coupling [29] of the colliding pair to a bound 
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singlet scattering length a s 
triplet scattering length a t 


52.61 a 
64.41 a 


vdW coefficient C% 
vdW length R vdw = i(2^C 6 /ft 2 ) 1/4 
vdW energy E vdw = h 2 /(2^,Rl dw ) 


2322 E h 
40.8 a 
h x 207.6 MHz 



Table 1. Important parameters of the interaction poten- 
tials of 6 Li- 40 K, taken from the potentials of Ref. [30]. The 
van der Waals (vdW) parameters describe the long range 
part of the potential, —Ce/r 6 , where r is the interparticle 
distance. Here, ao = 0.529177 x 10~ ln m is the Bohr radius 
and E h = 4.359744 x 10~ 18 J represents a hartree. 
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Fig. 2. (Color online) s-wave scattering properties of the 
ac channel, as a function of magnetic field. A coupled 
channels calculation (CC, solid line) is compared to a fit 
using Eqs. (2)-(4) (dashed line). The top panel shows the 
real part of the scattering length, while the two-body loss 
rate coefficient is shown in the lower panel. 



state of the same Mp = mi + m,2- Here mi 2 are the pro- 
jections along the magnetic field quantization axis of the 
total angular momenta of atoms 1 and 2, / 1,2, and Mp is 
that of the total spin angular momentum, F = fi + f%. 
We note that / and F are only good quantum numbers 
at zero magnetic field. 

For a pair of atoms in an excited Zeeman channel, there 
are two processes that can cause two-body collisional loss. 
Spin-exchange coupling can lead to inelastic spin relax- 
ation (ISR), in which the colliding pair is coupled into 
an energetically lower channel of the same Mp and the 
same partial wave £. Since each resonance considered in 
this work is in the energetically lowest s— wave channel 
of the relevant Mp, ISR does not occur. Spin dipole cou- 
pling [29], however, can couple a colliding pair to chan- 



nels of a different Mp or I, under the constraints that 
Mtot = Mp + mi is conserved, and the change in partial 
wave is given by At = ±2. Here, mg is the projection of 
& along the magnetic field quantization axis. For the res- 
onances considered here, spin dipole coupling links s- and 
d- waves, d- and g- waves, etc., with odd partial waves ex- 
cluded by symmetry requirements. For the ac channel, the 
two main decay pathways are the aa and ab d— wave chan- 
nels. Consequently, a basis including all s— and 0?— wave 
channels with M tot = —2 is sufficient. 

Scattering properties in the vicinity of the resonance 
at 154.75 G are shown in Fig. 2, along with the fit from 
Eqs. (2)-(4). The calculation used a collision energy of 
&b x 1 pK, while the fit assumes zero temperature. The fit 
gives excellent agreement, with only small deviations visi- 
ble outside the core of the resonance, and at the very cen- 
ter, where effects related to finite collision energies become 
important. The background scattering length near the res- 
onance is 63.0 ao, a suitable value for evaporative cooling, 
while the resonance width of 0.88 G makes it easily ac- 
cessible experimentally. Suppression of collisional losses is 
provided by the /cb x 14 mK height of the d-wave bar- 
rier being greater than the Zeeman splitting (fee x 1.8 mK 
for ab, &b x 3.6 mK for aa) between the entrance and 
exit channels. Consequently, decaying pairs must tunnel 
through the centrifugal barrier. The resulting resonance 
length is 4.0 x 10 6 ao- This is comparable to results we have 
found for much broader, entrance-channel dominated res- 
onances, such as the ee resonance of 85 Rb at £?o = 155 G, 
which has 10 6 a . We note that three-body 

effects, not included in our calculations, arc also of signif- 
icance for experiments. 



3.2 Scattering properties: Experiment 

3.2.1 Experimental conditions 

The basic procedures to prepare the Fermi-Fermi mix- 
ture near the 155 G Feshbach resonance are described in 
Ref. [8]. Here we briefly summarize the main experimen- 
tal parameters, and mention some issues of particular rel- 
evance for the present experiments. 

Our optical dipole trapping scheme employs two stages. 
In the first stage, we use a high-power laser source (200 W 
fiber laser) to load and evaporatively cool the mixture 
[6, 8] . As the quality of this high-power beam suffers from 
thermally induced effects such as spatial shifts and ther- 
mal lensing effects, we transfer the mixture into a second 
trapping beam that uses less laser power and is optimized 
for beam quality. This beam serves as the trapping beam 
in the second stage where all the measurements are per- 
formed. As the laser source we either use a broadband 
5W fiber laser (IPG YLD-5-1064-LP, central wavelength 
1065 nm) or a 25 W single-mode laser (ELS VersaDisk 
1030-50, central wavelength 1030 nm) 2 . In both cases the 



2 Specific product citations are for the purpose of clarifica- 
tion only, and are not an endorsement by the authors, JQI or 
NIST. 
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trapping potential (Gaussian beam waist 41 fim) is essen- 
tially the same, but we found that the broadband fiber 
laser can induce inelastic losses [39] . For our measurements 
on elastic interactions (Sec. 3.2.2) we used the broadband 
laser, and we later switched to the single-mode laser for 
the measurements of inelastic decay (Sec. 3.2.3). At a laser 
power of 70 mW the trapping frequencies for Li (K) are 
13 Hz (4.5 Hz) axially [40] and 365 Hz (210 Hz) radially, 
and the trap depth is 1.6 /iK (3.6 /iK). 

The mixture contains about 10 5 Li atoms at a tem- 
perature T Li w 140 nK together with about 2 x 10 4 K 
atoms at a temperature T K rj 160 nK; the temperatures 
are measured by time-of-flight imaging. Note that the two 
species are not fully thermalized at this point, such that 
T K > T Ll . In terms of the corresponding Fermi tempera- 
tures Tp l = 490 nK and Tf = 140 nK, the temperatures 
can be expressed as T^/Tp « 0.3 and T K /Tf « 1.1. 

The magnetic field was calibrated by driving RF tran- 
sitions between the b and the c state of K and using the 
Breit-Rabi formula [41]. In our set of measurements on 
the elastic scattering properties (Sec. 3.2.2) the magnetic- 
field uncertainty was about 20 mG, with a substantial con- 
tribution from magnetic field ripples connected with the 
50-Hz power line. In the later experiments on inelastic de- 
cay (Sec. 3.2.3) we could reduce this uncertainty down to 
about 5 mG. 



3.2.2 Elastic Scattering 

Our measurements on elastic scattering are based on the 
observation of sloshing motion, serving as a simple and 
sensitive probe for interspecies interactions [42-45] . With- 
out interaction both components would oscillate indepen- 
dently with their different sloshing frequencies. The inter- 
action induces friction between the two components and 
thus leads to damping. In the regime of weak interactions 
with up to a few scattering events per oscillation period, 
the damping rate can be assumed to be proportional to 
the elastic scattering cross section. Note that an alterna- 
tive approach, based on cross-dimensional relaxation, was 
followed in Ref. [9]. 

Here we restrict our attention to the slow axial sloshing 
motion. We excite this motion by an additional infrared 
beam intersecting our trapping beam [46]. The magnetic 
field is quickly ramped to the final setting that is applied in 
the measurements. After a variable hold time in the trap, 
we image both clouds to record their damped oscillatory 
motions. Our data analysis is based on the position of the 
K cloud, which is completely immersed in the much larger 
Li cloud. Its motion is analyzed by fitting a simple damped 
harmonic oscillation 



z(t) = Ae c * sin (ut + 4>) + z 



(8) 



to the observed axial center-of-mass position. Here A is the 
oscillation amplitude, C represents the damping rate, uj is 
the oscillation frequency, and zq the equilibrium position. 

Near the Feshbach resonance, the observed damping 
strongly depends on the magnetic field, as demonstrated 
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Fig. 3. (Color online) Samples of the damped sloshing 
motion of 40 K for three different settings of the magnetic 
field (upper panel 154.60 G, middle panel 154.96 G, lower 
panel 155.66 G). The solid circles represent the experimen- 
tal data, with uncertainties smaller than the size of the 
symbols. The solid lines are fits under the assumption of 
a simple damped harmonic oscillation. 



by the three sample oscillations displayed in Fig. 3. The 
measured damping rate as a function of the magnetic field, 
shown in Fig. 4, reflects the characteristic Fano profile of 
the elastic scattering cross section. The measured rates 
vary over three orders of magnitude, prominently showing 
both the pole of the resonance and its zero crossing. 

We analyze the observed magnetic-field dependence of 
the damping under the basic assumption that the rate 
C is proportional to the elastic scattering cross section, 
which itself is proportional to a 2 . Moreover, we take a 
background damping [47] into account which is indepen- 
dent of the interspecies interaction and express the total 
magnetic-field dependent damping rate as 



C(B) = Aa(Bf + Co 



(9) 



For the scattering length a(B) we use the result of the 
coupled-channels calculation as presented in Sec. 3.1. The 
theory has an uncertainty in the resonance position of 
the order of 100 mG, limited by the accuracy of the spec- 
troscopically derived potentials. We therefore allow for a 
magnetic-field offset by setting 



a(B) = a cc (B + S) 



(10) 



where a cc (B) refers to the scattering length resulting from 
the coupled-channels calculation (Sec. 3.1) and S is used 
as a free parameter. Based on Eqs. (9) and (10) we fit the 
experimental damping data with the three free parameters 
A, Co, and 8. 

The fit result, shown by the solid line in Fig. 4, shows 
excellent agreement with the experimental data. For the 
background damping of the non-interacting mixture, the 
fit yields Co = 0.0053(3) s _1 , which is consistent with inde- 
pendent measurements on K without Li [47] . For the mag- 
netic field offset parameter, the fit yields S = +69(3) mG. 
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Fig. 4. (Color online) Elastic scattering near the 155 G 
interspecies Feshbach resonance. The measured rates for 
the 40 K sloshing motion are shown (filled red circles) to- 
gether with a fit based on a coupled-channels calculation 
of the scattering length (solid black line) . For comparison 
we also plot damping rates measured for the non-resonant 
ab channel (solid blue squares) together with a theoretical 
line derived from the corresponding non-resonant scatter- 
ing length (dashed blue line). The error bars shown for the 
experimental data indicate the fit errors of the damping 
rate. 

Based on the theoretical results of Sec. 3.1 and this shift, 
we obtain a resonance position of 154.69(2) G with the 
20 mG calibration uncertainty being the dominant error 
source. 

The experimental data can also be analyzed indepen- 
dently of the coupled-channels calculations by using the 
standard Feshbach resonance expression (Eq. (3) in the 
limit 7b — > 0) for a fit in which the width A is kept as a 
free parameter. Our corresponding result A = 920(50) mG 
is consistent with the prediction A = 880 mG resulting 
from the coupled-channels calculation. 

For comparison, we have also investigated elastic scat- 
tering in a Li-K mixture in the ab channel (solid squares in 
Fig. 4). which near 155 G is weakly interacting. Our mea- 
surements show a damping of the sloshing motion that 
is consistent with the predicted non-resonant scattering 
length of 68 ao for this channel (solid line) . 

3.2.3 Inelastic Scattering 

To probe inelastic decay, we first prepare a weakly in- 
teracting, long-lived Li-K mixture in the ab channel at a 
variable magnetic field near 155 G. We then apply a short 
RF 7r-pulse (duration 90 /xs) to quickly transfer the mix- 
ture into the ac channel. This transfer method avoids fast 
magnetic field ramps and thus any waiting time for the 
magnetic field to be settled to a constant value before 
measurements can be taken. 
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Fig. 5. (Color online) Decay of K immersed in the Li 
cloud for different magnetic fields, far away from reso- 
nance at 154.766 G (triangles), very close to resonance 
at 154.712 G (circles), and for a setting in between at 
154.727 G (squares). The scatter of the data indicates the 
shot-to-shot variations of the measurements. As the life- 
times are plotted on a logarithmic scale, the good linear 
fits indicate pure exponential decay. 



Figure 5 shows example decay curves. The K loss is 
essentially exponential, which results from the fact that 
the K cloud is immersed in a much larger Li sample [6]. 
In this regime, the Li cloud serves as a large bath with 
essentially constant density. Here the loss curves do not 
allow us to distinguish between two-body losses where one 
K atom interacts with one Li atom and such three-body 
losses, where one K atom interacts with two Li atoms. 
Three-body losses resulting from two K atoms interacting 
with one Li atom would not lead to exponential loss. 

We analyze the decay under the hypothesis of dom- 
inant two-body loss, which is motivated by the decay- 
ing character of the Feshbach resonance as discussed in 
Sec. 3.1. The total K decay rate can be approximated by 
r = Kzfou) + i~bg, where fbg is a small background loss 
rate [48]. The mean Li number density is given by (nu), 
where the angle brackets denote averages weighted by the 
K density distribution. For our experimental parameters 
we obtain (n-Li) = 5.9 x 10 11 cm -3 , which is about 75 % of 
the peak density in the center of the Li cloud. 

Figure 6 shows the measured values for the loss rate co- 
efficient K2- The data show the expected resonance behav- 
ior (Sec. 3.1). The values peak at the center of the Fesh- 
bach resonance and strongly decrease within a few 10 mG 
away from the center. For large scattering lengths, the 
data follow the scaling oc a 2 according to Eq. (7). The 
observed resonance behavior thus confirms our assump- 
tion of the dominant two-body nature of loss. Three-body 
losses in a two-component Fermi mixture would show a 
much stronger dependence on a [49], not consistent with 
this observed behavior. However, significant three-body 
loss contributions may be present very near to the reso- 
nance center. 

The figure also shows three theoretical curves, cal- 
culated for three different values of the collision energy 
(&b x 1 pKj representing the zero energy limit, fee x 100 nK, 
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Fig. 6. (Color online) Inelastic loss near the 155 G in- 
terspecies Feshbach resonance. The measured values for 
the two-body loss rate coefficient K 2 (solid circles) are 
compared with theory. The three theoretical curves (solid 
lines) represent three different collision energies E/k-Q 
(1 pK, labelled to indicate the zero energy limit, 100 nK, 
and 300 nK), showing the limiting effect of finite collision 
energy. The error bars represent the statistical errors from 
fitting the loss curves. 



and fee x 300 nK) in a range relevant for our experiments. 
As a typical value for the collision energy, we can consider 
an estimate of 200 nK [50] . Very close to the resonance the 
theory curves illustrate how K 2 increases in the zero tem- 
perature limit up to a value corresponding to b = 2a res - In 
the case of non-zero collision energies it is limited to lower 
values. For magnetic detunings exceeding about 20 mG, 
the effect of the finite collision energies can be neglected in 
the interpretation of the experimental data, which makes 
the comparison between theory and experiment straight- 
forward. Here we find excellent quantitative agreement, 
confirming two-body decay as the dominant loss mecha- 
nism. Very close to the center of the resonance the situa- 
tion is more complicated. If one completely attributes loss 
to two-body decay, the 100 nK curve provides an excellent 
fit to the experimental data. This, however, is somewhat 
below our estimate of 200 nK for an effective collision en- 
ergy, which may point to additional three-body losses at 
the very center of the resonance. 

To extract the precise resonance position we proceed 
in an analogous way as for analyzing the elastic scattering 
data, allowing for a small magnetic field shift 8 between 
theory and experiment. We write the actual loss coeffi- 
cient as K 2 (B) = K 2: cc{B + 8), where K 2: cc refers to the 
coupled-channels result for K 2 as discussed in Sec. 3.1. 
In the fit, we exclude the three experimental data points 
that exceed 3 x 10~ u cm 3 /s to avoid the region where fi- 
nite collision energies become important. This also makes 




Fig. 7. (Color online) Molecular binding energy as a 
function of magnetic field. MQDT (blue line) refers to 
the three-parameter model of Rcf. [51], while CC (green 
points) indicates a coupled channels calculation. Suffi- 
ciently close to the resonance the binding energy converges 
to the universal result of Eq. (11) (red dashed line). 



sure that the loss data are dominated by two-body decay. 
The shift 8 is the only free parameter, and we obtain a 
small value of 8 = +38(1) mG, well in the range of the 
theoretical uncertainty. We finally obtain a resonance po- 
sition of Bq = 154.707(5) G, where the main uncertainty 
results from the magnetic field calibration. Within the ex- 
perimental uncertainties this value is consistent with the 
less precise resonance position obtained from elastic scat- 
tering measurements. 



3.3 Bound state properties 

In the context of Feshbach molecules, universality refers to 
the range of magnetic fields sufficiently close to resonance 
within which the molecular and scattering properties can 
be described solely by the atomic masses and the scatter- 
ing length a(B). Within this region, the molecule has the 
form of a halo state, in which a significant part of the wave- 
function lies beyond the classically allowed outer turning 
point of the potential. This results in a strong enhance- 
ment of the lifetime of a decaying bound state [36,37]. The 
universal binding energy is given by 



Et 



2^a(B) 2 



(11) 



Calculations of the relevant bound state energy using the 
coupled channels method and the simplified three-parameter 
model of Ref. [51] are shown in Fig. 7. We note that the 
three-parameter model, while useful for bound state and 
resonance characterisation, does not couple partial waves 
and so can not be used for calculating decay properties in 
the present case. The energy variation is linear for bind- 
ing energies greater than a few tens of kHz, having a rela- 
tive magnetic moment of 8/j,/h = 2.3MHz/G with respect 
to the ac threshold. The universal region, as can be seen 
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Fig. 8. (Color online) Calculated molecular lifetime as a 
function of magnetic field. The molecular lifetime varies 
at lower field according to its tunnelling rate through the 
<i-wave centrifugal barrier into the exit channels. The ar- 
rows indicate the sharp increase in lifetime as the univer- 
sal regime is entered, and the region away from resonance 
where the molecular lifetime converges to that of the bare 
resonance state. 

from the inset of Fig. 7, covers a magnetic field range of 
order mG. This makes it hard to access experimentally. 
The universal region is wider for broad, entrance channel 
dominated resonances [29]. However, in the present case, 
the suppression of decay by the centrifugal barrier allows 
the molecules to have a long lifetime in the nonuniversal 
regime. 

We now consider the lifetime of 6 Li- 40 K molecules close 
to the ac resonance at 155 G. Outside the very narrow uni- 
versal region, the analytic approach of Ref. [37] does not 
apply. We therefore derive the molecular lifetime from a 
coupled channels scattering calculation including the two 
open d-wave channels into which it decays. The spin-dipole 
induced decay discussed in the previous section is medi- 
ated by the bound state causing the resonance. Around 
the bound state energy, the off-diagonal T matrix element 
between the two decay channels follows the form 

|T 12 | 2 = ^ — g . (12) 

(E - E Ies f + h 2 (2i±2i) 

Here, 71 /2 and 72 /2 are the decay rate of the molecule into 
the ab and aa channels, respectively. The total decay rate 
is given by (71 + 7a)/2. This calculation includes the en- 
trance channel component and so reproduces the increase 
in lifetime as the Feshbach molecule takes on a halo form. 
This should be distinguished from the decay rate of the 
bare resonance state which appears in Eqs. (3) and (4). 

Our calculated lifetimes are shown in Fig. 8. As dis- 
cussed above, a sharp increase in lifetime occurs as the 
universal region near Bq is approached. Above the max- 
imum lifetime shown in the Fig. 8, the decay peak de- 
scribed by Eq. (12) narrows to the point where we can 
no longer resolve it in our calculations. The slower in- 
crease as B is moved away from Bq occurs because the 




B(G) 

Fig. 9. (Color online) Energies of the bound states under- 
lying the resonances studied in this paper, as a function 
of magnetic field. Collision thresholds of relevant chan- 
nels are shown as dotted lines, labelled at the right of the 
figure. The resonances arise from three zero-field bound 
states of F = 4, 5 and 6 (labelled at the left of the fig- 
ure), which are projected into channels of |M to t| < F at 
non-zero field. Deeper bound states are not shown for rea- 
sons of clarity. Avoided crossings between bound states of 
the same M to t give rise to three identifiable groups of res- 
onances, indicated by the same symbols used in Fig. 10 
and Table 2. Within each group, resonance parameters 
vary smoothly (see Fig. 10 and text). 

bound state moves further behind the centrifugal barrier. 
Decay from tunnelling through the barrier is then fur- 
ther suppressed. The lifetime of molecules in the vicinity 
of the 155 G resonance was measured by Voigt et al. [5]. 
They observed a sharp increase in lifetime near the res- 
onance, with which our results qualitatively agree. Their 
measured background lifetime of ~ 3 ms away from res- 
onance is lower than our calculated minimum of 11ms. 
However, our calculations do not include relevant atom- 
dimer and dimer-dimer collisions, and so may be consid- 
ered as an upper bound to experimentally observable life- 
times. A lifetime of several ms will permit measurements 
and manipulation of the Feshbach molecules. 

4 Survey of resonances 

In this section we discuss resonances occurring in various 
channels of the 6 Li- 40 K mixture. We focus on channels 
with 6 Li in the a state and 40 K in the lower (/ = 9/2) man- 
ifold, for which inelastic spin-exchange collisions do not 
occur. At zero magnetic field, there are three bound states 
of F = 4, 5 and 6 in the range 200 MHz to 300 MHz below 
these thresholds, as shown in Fig. 9. At nonzero magnetic 
field, these states are projected into their Zeeman sub- 
levels, which give rise to Feshbach resonances when degen- 
erate with the collision threshold of a channel of the same 
M to t- Consequently, three proximate resonances are found 
in channels of —3 < M to t < 4. The bound state underly- 
ing each resonance adiabatically correlates with one of the 
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Fig. 10. (Color online) Theoretical survey of 6 Li- 40 K s- 
wave resonances. Each panel shows a resonance parameter 
as a function of M tot : the width A (top), strength s res 
(middle), and length a rcs (bottom). The 6 Li atom is in 
the a state, except for M to t — —5 for which it is in b. 
The 40 K atom is in the lowest Zeeman state producing 
the relevant M tot . The symbols, also used in Fig. 9 and 
Table 2, correspond to the resonance groups discussed in 
the text. 



zero field F states. We note that the bound state energies 
shown in Fig. 9 were produced with the three-parameter 
model of Ref. [51], which produces slightly different reso- 
nance locations to the coupled channels calculations that 
follow. 

We have performed a coupled channels analysis of each 
resonance, analogous to that performed with the asymp- 
totic bound state model in Ref. [7]. With our more rigor- 
ous approach [31], we obtain good agreement with all ex- 
perimental data, including the new set of measurements on 
the 155 G resonance presented in Sec. 3.2. The simplified 
approach of Ref. [7] seems to underestimate the widths 
of the resonances by almost a factor of two. With the 
coupled channels approach, we can also study the decay 
properties of the resonances. The resonance parameters 
are shown in Fig. 10, and tabulated in Table 2. We group 
the resonances of each channel that are lowest (O)i mid- 
dle (□) and highest (0) in Bq, using the indicated symbols 



to distinguish the resonance groups in Figs. 9 and 10, and 
Table 2. Within each of these groups, resonance proper- 
ties vary smoothly as a function of M to t- The resonances 
with Mtot = — 4 and 5 have properties consistent with 
the lowest and highest group, while the ba resonance with 
Mtot = —5 has substantially different properties. This is 
due to F being a good quantum number only at zero mag- 
netic field, and several bound states having avoided cross- 
ings in the relevant range of magnetic field. There are 
several resonances with A > 1 G, offering good opportu- 
nities for control of collisional properties. However, several 
other factors are also useful for deciding the suitability of 
a resonance for a given application. 

One parameter used for quantifying the extent to which 
a resonance is entrance-channel dominated is the reso- 
nance strength parameter [29], defined by 



dbg SfiA 
1l W 



(13) 



Here, a — [4-7r/r(l/4) 2 ]i? vd w ~ 0.956i? v dw is the mean 
scattering length [52], and E = h 2 /(2^a 2 ) « 1.094Evdw 
is the associated energy If ,s rcs > 1, the bound state and 
near-threshold scattering states are concentrated in the 
entrance channel over a magnetic field range comparable 
to A. In the present case, all resonances are closed-channel 
dominated, as shown in the middle panel of Fig. 10. The 
background scattering lengths of the resonances are all in 
the range 55 ao to 70 ao, and the relative magnetic mo- 
ments are in the range 2MHz/G to 2.6MHz/G. Conse- 
quently, the resonance strength follows trends similar to 
the resonance widths. For the ac 155 G resonance we have 
s r es = 0.014. This is reflected in the universal region being 
only a few mG wide, as discussed above. 

Calculated resonance lengths are shown in the lower 
panel of Fig. 10. The range in a ros is approximately 3 
orders of magnitude, with better stability in channels of 
lower M to t- This occurs because the energy gaps between 
higher channels arc larger, reducing the height of the cen- 
trifugal barrier through which the decaying atoms tunnel. 
However, all the resonance lengths are sufficiently high 
that we expect each resonance with A > 1 G to be very 
useful. 

In view of interaction control in a strongly interacting 
gas, we now discuss three selected resonances that have 
received particular attention in experiments: the 168 G 
resonance in the aa channel [8], the 155 G resonance in 
the ac channel [5] (see Sec. 3), and the 114 G resonance 
in the aj channel. In practice, the possible degree of con- 
trol is limited by uncertainties (drifts and fluctuations) of 
the magnetic field. A corresponding figure of merit is the 
maximum controllable scattering length a c tri = ^bgA/Ss, 
where 5b stands for the magnetic field uncertainty. As- 
suming a realistic value of 5b = 5mG one obtains o c t r ] ~ 
1600 ao, ll'OOOao, and 21'OOOao for the three resonances 
considered (168 G, 155 G, and 114 G, respectively). On one 
hand, this can be compared with the typical requirement 
°f l a l ^ 5000 ao for attaining strongly interacting condi- 
tions. On the other hand, it can be compared with the con- 
dition for universal behavior |a| 3> a/s ros [53], which re- 
quires a ctr i ^> 17'OOOao, 2800 ao, and 1450 ao- This shows 
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Table 2. Survey of s-wave resonances in 6 Li- 40 K. The first two columns give the channel a/3 and total angular 
momentum projection M to t, with a and (3 representing the Zceman state of Li and K, respectively. The third column 
gives the symbol for the corresponding group of resonances that is used in Figs. 9 and 10. The next three columns 
give experimental values of the resonance location Bq and width A, where available, with references. The remaining 
columns give the results of the coupled channels calculations performed for the current work - Bq and A, as well as 
the background scattering length ah g , relative magnetic moment 5[i, resonance length o res , resonance strength s rcs , 
and decay rate in magnetic field units, 75. Note that a rcs is not defined for the stable aa channel. Note also that the 
experimental values for Bq from Ref. [4] are subject to typical uncertainties of about 0.5 G. 



that the resonance at 168 G is too narrow for control- 
ling a strongly interacting Fermi-Fermi mixture, but the 
other resonances at 155 G and 114 G are broad enough. 
Although the 114 G resonance is wider than the 155 G res- 
onance by a factor of 2.1, inelastic loss is 3.7 times faster. 
The higher collisional stability is an important advantage 
of the 155 G resonance. 



5 Conclusions 

We have characterized elastic and inelastic scattering near 
Feshbach resonances in the 6 Li- 40 K mixtures. The pres- 
ence of open decay channels for all broader resonances has 



two important consequences. Atomic two-body collisions 
acquire a resonantly enhanced inelastic component, which 
unavoidably limits the stability of an atomic Fermi-Fermi 
mixture with resonantly tuned interactions. When Fesh- 
bach molecules are created via these decaying resonances, 
they will undergo spontaneous dissociation. 

The intrinsic decay has important consequences for 
present experiments towards strongly interacting Fermi- 
Fermi mixtures. Under typical experimental conditions, 
the lifetime of a Fermi- Fermi mixture with resonantly tuned 
interactions (a — » ±00) will be limited to ~10 ms. This 
in general means a limitation of possible experiments to 
short time scales, such as the observation of the expansion 
of the mixture after trap release [54, 55] or measurements 
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of fast collective oscillation modes [56-58]. Experiments 
that require long time scales, such as precise studies of 
equilibrium states [59,60], may be problematic in this de- 
caying mixture. 

The short lifetime of the Feshbach molecules, also be- 
ing of the order of 10 ms, excludes the production of a long- 
lived molecular Bose-Einstein condensate (mBEC) such as 
formed in 6 Li [61,62]. Transient ways to form mBECs, as 
demonstrated for 40 K [63], will still be possible. The de- 
tection of fermionic condensates by rapid conversion of 
many-body pairs into molecules [64] also seems to be a re- 
alistic possibility. Moreover, the predicted increase of the 
molecular lifetime for larger binding energies can be of 
general interest for the coherent manipulation of Feshbach 
molecules [65] and in particular for optimizing the start- 
ing conditions for a transfer to the ro- vibrational ground 
state [66-68]. 

Finally, the question of which Feshbach resonance pro- 
vides optimum conditions for interaction tuning in the 
6 Li- 40 K has no straightforward answer. All the broad res- 
onances occurring in the channels ac-aj (widths 0.88 G - 
1.97 G) seem to be well suited for controlled interaction 
tuning. Because of the tradeoff between width and stabil- 
ity, the best choice will depend on the particular applica- 
tion. 
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and the European Science Foundation (ESF) within the Euro- 
QUAM/FerMix project and support by the FWF through the 
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